Rapid non-adiabatic loading in an optical lattice 
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We present a scheme for non-adiabatically loading a Bose-Einstein condensate into the ground 
state of a one dimensional optical lattice within a few tens of microseconds typically, i.e. in less 
than half the Talbot period. This technique of coherent control is based on sequences of pulsed 
perturbations and experimental results demonstrate its feasibility and effectiveness. As the loading 
process is much shorter than the traditional adiabatic loading timescale, this method may find many 
applications. 
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Numerous works related to optical lattice trapping 
have been published, especially using coherent atoms or 
Bose-Einstein condensate (BEC) [3|, since it has vari- 
ous applications in quantum computation, simulation of 
basic condensed matter physics, atomic clocks, etc. A 
common concern to those experiments is how to load a 
BEC into the lattice without excitation or heating Q-Q. 
In most cases, one chooses to turn on the light field adia- 
batically and the loading process usually lasts up to tens 
of milliseconds, during which one tries to avoid excita- 
tions to higher states. A long loading may be problem- 
atic for quantum computation experiments as it increases 
the time during which decoherence can occur and reduces 
speed when atoms stored in optical lattices shall be in- 
terrogated and reloaded several times. Furthermore, adi- 
abatic loading is difficult to obtain near the critical point 
of phase transitions for finite temperatures even well bel- 
low the BEC critical temperature. 

An alternative idea is to prepare the BEC in the ground 
state of the lattice, and then turn on the lattice suddenly. 
This 'preparing' process can be much shorter compared 
to the adiabatic loading (about 30 /is, as we can see 
later in this article). One way for realizing this 'prepar- 
ing' process can stem from nonholonomic coherent con- 
trol [l|, in, 0- In this approach, a sequence of two well 
chosen Hamiltonians are imposed on the system and the 
preset duration of each step is modified according to a 
defined cost function, in order to get the aimed evolution 
operator as well as the target state. 

Our proposal for designing and computing the pre- 
loading process is reminiscent of this method, but we fo- 
cus more on the feasibility of its experimental implemen- 
tation, rather than the stringency or universality of the 
approach. Experimentally, we achieve the non-adiabatic 
loading with a very few pulses using directly the results of 
our computational design. According to our analysis, this 
method would be valid on many other occasions, such as 



loading the condensate directly into excited states of the 
lattice. We first introduce a general method for deter- 
mining the sequence of steps to be applied to the system 
before giving an account of our experimental results. 

Suppose that before the sudden loading at time to of 
the lattice with optical depth Vb, m steps have been ap- 
plied. The ith step corresponds to a Hamiltonian Hi kept 
constant during ti. The final state |'0(to)) is given by: 



mo)) = n ^'i^o) 
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where Ui = e~*^^**/^ is the evolution operator of the ith 
process. With l^pa) the aimed state, which is the ground 
state of the lattice in the context of this paper, the total 
population of the excited states A^e at time to, is: 



= 1 - \{xPamto))\' 
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Our goal is to properly choose Hi and ti so that A^e is 
small enough to be neglected in the experiment, or the 
final state |'0(to)) to be close enough to 

The target state \tlja) can be obtained simply if atomic 
interactions are neglected. The effect of the optical lat- 
tice on the atoms corresponds to a stationary periodi- 
cal potential with periodicity Al/2, spontaneous emis- 
sion being neglected since the laser used for the lattice 
is sufficiently far detuned. So neglecting atomic interac- 
tions, the Hamiltonian of the interaction of atoms with 
the optical lattice is: 
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where m is the atom mass, kj, is the wave number of lat- 
tice light and Vq the lattice depth and p is the momen- 
tum of each atom. The corresponding Talbot time [8] is 

Tt = ^j^, where ujr = is the one photon recoil fre- 
quency. Because of the periodicity of the potential, {ipa) 
can be decomposed over a reduced basis of plane waves 
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\2nhkL). One obvious choice for each Hi is to take the 
Hamiltonian corresponding to the interaction of atoms 
with a standing wave with same periodicity Al/2 so that 
the Hilbert space is hmited to the basis of states \2nhkL)' 
For this purpose, the power of the same laser as the one 
used for the final lattice loading is simply adjusted and 
each Hamiltonian Hi is obtained from Eq. ([3j) after sub- 
stitution of V{) by the new lattice depth Vi. 

More precisely, as Hi has spatial periodicity, we 
get its eigenstates by solving the equation Hi\n^q) = 
En.qWiO) [9] 7 where \n^q) is a Bloch state with n the 
band index and q the quasi momentum and E^^q is 
the corresponding eigenenergy. We use the notation 
\n^q^Vi) for denoting the Bloch states for a Vi lattice 
depth. Since only Bloch State with = is pop- 
ulated initially, no other quasimomenta can be popu- 
lated during the sequence of pulses. Thus the state of 
the system can be spanned over the momentum eigen- 
states basis \2ihkL)^ independent on the potential depth 
Vi and the evolution operator can be written as the fol- 
lowing matrix: lJi{Vi,ti) = C{Vi)E{Vi,ti)C{Viy , where 
C{Vi) is the unitary matrix of transition between the 
Bloch States basis and the momentum eigenstates ba- 
sis with matrix elements G{Vi)£ri = {2ihkL\n^q = O^Vi) 
and where E(T^,t^) is a diagonal matrix with elements 
'E,{Vi,ti)nn = exp{-iEn,q=o{Vi)ti/h). Bccausc of the 
simple form of the potential, these matrices are easily 
obtained and from those the wave function's evolution is 
solved numerically. 

A traditional simple approximation to solve this prob- 
lem of a standing wave of light pulse interaction with 
atoms is to omit the effect of kinetic energy. As shown in 
the works [H, [lo| , the motional term in the Hamiltonian 
may be neglected. For a square pulse, when the depth 
of the standing wave potential is V and the duration 
of the pulse is t, the wave function after the pulse can 
be expressed as ilj{t) = ?/;(0)e~*^ and decom- 

posed over different momentum states \2nhkL)^ using the 
Bessel functions. This approximation is valid only in the 
Raman-Nath regime, when the pulse is short enough that 
the displacement of the atoms is much smaller than the 
lattice period. From references [HI, [III, we can see exper- 
imental data deviate distinctly from the Bessel functions 
when the pulse duration is longer. For multi pulse cases, 
the motion term may be neglected during each pulse but 
taken into account during the intervals between pulses. 
This approximation is valid only when the sum of all 
thepulses' durations is far less than the Talbot time T^. 
In [8[, Tt = 10 /is, and two 100 ns pulses were applied, 
while the interval time can be up to 10 /is. However, 
as the total pulse duration 200 ns is far less than Tt, 
the whole process is still in Raman-Nath regime and we 
can see the theoretical prediction agrees well with the 
experimental results. 

We compare in Fig. 1 the results of the two meth- 
ods with our experimental data for the case of a single 
square pulse applied to the system. In our experiment, 
we first prepare a BEG of about 2 x 10^ ^^Rb atoms 
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FIG. 1: (color online) Relative atom number No in p — OhkL 
mode versus the pulse duration t of a single standing-wave 
pulse applied to the system. The black points with error 
bars are experimental data and the red dashed-dotted line is 
calculated with the simple method where the kinetic energy 
is neglected while the blue dashed line is given by the general 
method (see text). 



in the \F = 2,mi? = 2) hyperfine ground state in the 
magnetic trap, with longitudinal length L = 100 fim and 
width / = 10 /im [13]. The BEG is then loaded into 
a one-dimensional optical lattice with light wavelength 
Xl = 852 nm, along its axial direction. We suddenly 
turn on the standing-wave potential to a specific depth, 
hold the lattice for different durations release the con- 
densate and measure the relative atom number Nq in 
p = OhkL mode. Without considering the atomic motion, 
the Bessel function (red dashed-dotted line) failed to fit 
the experiment data when t is larger than about 10 /is. 
On the other hand, the numerical solution taking the mo- 
tional term into account (blue dashed line) predicts the 
experiment accurately up to 60 /is, although some other 
effects have been neglected, such as interaction between 
atoms and non-zero momentum width of the condensate. 
As a result, we will use the general method as it will im- 
pose less restriction on our sequence design. Note that 
Tt ^ 79 /is in our configuration and the fitted potential 
depth is IS Eji^ where Er = (hkLY /2m is the one photon 
recoil energy. 

In a first attempt to design the time sequence, we use 
four steps, each corresponding to the application of a 
potential with depth Vi and a duration ti as shown in 
Fig. 2a. Thus, from these eight free parameters Vi^ti 
with z = 1 to 4, we can obtain the evolution operators Ui 
and the exciting rate A^e according to our previous anal- 
ysis. If we want to get the minimum of A/'g, it requires to 
test all the combinations of the eight free parameters and 
need too much computing time. However, we find con- 
vergence to a local minimum is easy to obtain. For this 
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TABLE I: Time duration of each step ti (/xs) for different 
lattice depths Vo with unit Er. 
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FIG. 2: (color online) Two schemes for the non-adiabatic 
loading, (a) Four steps pre-loading sequence with both the 
potential depth and the duration of each step set as free pa- 
rameters, (b) Four steps pre-loading sequence with fixed po- 
tential depth and variable duration in each step, (c) Corre- 
sponding relative population A^e in excited Bloch states as a 
function of lattice depth, with red circles for scheme (a) and 
blue points for scheme (b). 



purpose we modify the parameters step by step, with a 
small change of one parameter every step, in order to get 
a smaller N^. If no smaller result for can be obtained 
by changing any of these parameters, the minimization 
of A^e is stopped and the corresponding combination of 
parameters is the result of our design. Note that differ- 
ent initial values of the parameters may result in different 
local minima. Actually we try more than one group of 
initial parameters in order to get a better result, though 
the optimized pulse sequences obtained for different and 
often lower lattice depths are already good input values. 
As shown in Fig. 2c, the exciting rate remains lower 
than 0.1%. The whole pre-loading process lasts for about 
30 /is, which is beyond the Raman-Nath Regime but still 
within the limit of our method. 

Nevertheless, it is not easy for this scheme to be ex- 
perimentally implemented. As the lattice depth varies 
with time, we need a feed back loop to control the laser 
power to get a stable optical lattice as the atoms are very 
sensitive to fluctuation of the laser power. We also need 
sharp rising and falling edges to ensure each step dura- 
tion is sufficiently similar to the designed one. If 0.5 /is is 
set for the upper limit of each rising and falling edge, the 
bandwidth of the power controlling loop should be larger 
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than 1 MHz, which is hard to achieve. We would rather 
choose the scheme illustrated in Fig. 2b, where the step 
Vi and are fixed at Vq^ and the step V2 and V4 are 
fixed at 0, ti still being free parameters. This means no 
feed back loop but only a switch is needed so that the 
rising and falling time can be decreased to 200 ns eas- 
ily. In Fig. 2c, we can see the exciting rate is still lower 
than 0.1% for most lattice depths lower than 20 Er^ al- 
though some degrees of freedom of our sequence of pulses 
are locked. In TABLE I we show the designed time se- 
quences for several lattice depths. 

Furthermore, we can transfer the state \ipa) back to the 
initial state |?/^o) after releasing the lattice, using a similar 
process. TABLE I gives also the designed sequence, with 
time durations to and corresponding potentials, V5 
and being set to 0, while Vq and equal to Vb (see 
Fig. 3b). Note that this sequence is equivalent to the 
timely inverted pre-loading sequence if we ignore some 
minor deviations. 

To demonstrate the feasibility of our proposal, we did 
the experiment according to three different schemes. In 
the first scheme (Fig. 3a), we non-adiabatically load the 
BEC into the lattice, according to the designed process, 
hold it in the lattice for th = 10 ms and release the atoms. 
We can see the interference peaks, similar to the famil- 
iar pattern observed in adiabatic loading experiments, 
which indicates a successful loading without significant 
excitation and heating. In the second scheme (Fig. 3b), 
after the same loading and holding process, we use two 
additional pulses to transfer the atoms back to the origi- 
nal state IV^o)- Compared with the third scheme (Fig. 3c), 
where we directly turn on and off the lattice light without 
the pre-loading or post-releasing processes to prevent ex- 
citation but hold for the same period, the second scheme 
has little heating or disturbing effect on our BEC, which 
proves the effectiveness of our 'preparing' process of the 
lattice ground state. A small heating may still be ob- 
served at 12 Er. It may be due to interaction between 
atoms which are not included in our model. 

Fig. 4 shows the evolution of the wave function un- 
der a pulsed scheme shown in Fig. 2b. We can see that 
atoms can be transferred between different momentum 
orders only when the pulse is on. In the intervals, the 
amplitudes of different components keep constant, while 
the phases vary linearly. The phase of the :^AhkL order 
evolves four times faster than that of ^2hkL as the ki- 
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FIG. 3: (color online) Absorption images for three different 
time sequences. As the condensate was released from both 
the lattice potential and the magnetic trap, it expands freely 
for 30ms and then an absorption image is taken. For each 
time sequence, we show the results for two lattice depths of 
8Er and 12Er. 

netic energy is four times lager. After switching on the 
optical lattice {t = /is), the amplitude and phase of 
all components keep almost unchanged as the wave func- 
tion is in the ground state of the lattice. Note that the 
phase jump of the non-zero momentum components at 
t = — 27 //s is not significant since the amplitude of these 
components before that time is zero. 

By simply replacing the initial state Ii/jq) and the aimed 
state IV^a), our pulse sequence design is also applicable to 
other situations. For example, by setting \ipo) = \n = 
0, = 0, Vi) and {i/ja) = |n = 0, = 0, ¥2)^ we can get a 
process that can change the lattice depth from Vi to V2 
non-adiabatically without excitation. We can also load 
the BEG directly to symmetric excited states such as 
\n = 2^q = 0), by setting them as l^/^o). Our numerical 
simulation shows when Vq is around 12 Eji, there ex- 
ists theoretical possibility to load more than 99% of the 



atoms to |n = 2, g = 0). This state can be used for study- 
ing spontaneous transition from high lying Bloch bands 
to lower bands experimentally. If the optical lattice is 
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FIG. 4: (color online) Evolution of the different momentum 
components of the wave function. Amplitude (a) and phase 
(b) for the p = OhkL, =b2/i/cL and z\z4:hkL modes are shown 
with blue solid lines, red dashed lines and green dashed-dotted 
lines respectively. The black dotted lines in both figures in- 
dicate the corresponding time sequence of the standing wave 
potential. 



accelerated, anti-symmetric states and states with q 
could also be loaded. 

In general, our method can be used to produce states 
whose momentum components are discrete and separated 
equally by 2hkL' For example, we can divide one BEG 
into two momentum modes \ ^2mhkL) equally, with neg- 
ligible population in other orders. This technique is use- 
ful in atom interferometry and its theory has been de- 
veloped in reference [lo|, while neglecting the motional 
term in the Hamiltonian, as discussed earlier in this ar- 
ticle. Our method is not restricted by the Raman-Nath 
regime and thus gives more freedom for designing the 
pulse sequence. 

In conclusion, we proposed a method reminiscent of 
the nonholonomic coherent control technique for load- 
ing the BEG into one dimensional optical lattice non- 
adiabatically, within a much shorter loading time (usu- 
ally less than Tt/2) than the commonly applied adiabatic 
method (much longer than Tt). Our experimental results 
demonstrate the validity of this method. We claim that 
this numerical design process can be applied to various 
topics related to the interaction of a standing wave of 
light with ultra-cold gases. 
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